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I.   INTRODUCTION 

The  method  of  intermediate  problems  for  obtaining  lower  bounds 
to  eigenvalues  was  introduced  by  Weinstein  [l]1  as  a  means  of  esti- 
mating from  below  the  eigenf requencies  of  a  vibrating  clamped  plate. 
His  procedure  involves  constructing  a  solvable  base  problem  to  the 
given  problem  by  relaxing  boundary  conditions  so  as  to  remove  con- 
straints from  the  system.   The  eigenvalues  of  this  base  problem 
thus  give  lower  bounds  to  the  eigenvalues  of  the  given  problem. 
Solvable  problems  intermediate  between  the  base  problem  and  the 
given  problem  are  constructed  by  systematically  restoring  the  re- 
moved constraints.   The  eigenvalues  of  these  problems  give  lower 
bounds  improving  toward  the  eigenvalues  of  the  given  problem. 
The  method  has  subsequently  been  formulated  in  terms  of  the  theory 
of  operators  in  Hilbert  spaces  [2,3]. 

A  different  procedure  for  constructing  intermediate  problems 
has  been  given  by  Aronszajn  [d].   In  this  procedure,  a  base  oper- 
ator, whose  eigenvalue  equation  can  be  solved,  is  constructed 
from  a  given  operator  by  subtracting  a  positive  definite  part. 
The  eigenvalues  to  the  base  operator  are  lower  bounds  to  those  of  * 
the  given  operator.   By  means  of  a  projection  technique  operators 


Numbers  in  brackets  refer  to  the  literature  references  located 
after  the  appendices. 


smaller  than  the  subtracted  operator  are  added  to  the  base  operator 
to  give  solvable  intermediate  operators  whose  eigenvalues  approach 
those  of  the  given  operator  from  below. 

Recently,  Baz.ey  and  Fox  [5l  6,  7]  have  used  the  projection 
technique  of  Aronszajn  to  construct  intermediate  operators  which 
may  be  solved  or  approximated  from  below  in  terms  of  finite  matrix 
eigenvalue  problems.   Their  methods  are  applicable  to  many  of  the 
Hamiltonians  of  atomic  and  molecular  systems. 

In  the  first  part  of  this  dissertation  we  develop  the  procedure 
of  Baz.ey  and  Fox  for  constructing  intermediate  Hamiltonians  along 
with  three  of  their  methods  of  solution.   We  also  introduce  another 
method  of  solution  which,  while  possessing  undesirable  features, 
appears  to  offer  some  advantages  in  determining  lower  bounds  to  the 
ground  states  of  atomic  and  molecular  systems. 

•n  the  last  part  we  apply  these  procedures  to  the  helium  atom 
and  to  a  guasi-atomic  system  consisting  of  two  electrons  confined 
within  an  isotropic  potential  well. 


II.   INTERMEDIATE  HAMILTONIANS2 


We  consider  a  time  independent  Schrt/d  i  nger  equation, 

gt4  =  £  <J     ,  (2-0 

for  which  the  Hamiltonian  at     can  be  decomposed  into  two  parts, 

M  =   K°  +  V  ,  (2.2) 

such  that  V  is  positive  definite,  and  ^*t  has  known  eigenvalues 
and  e igenf unct ions ,  is  bounded  from  below,  and  possesses  a  spectrum 
consisting  of  at  least  an  initially  discrete  set  of  levels. 

Since  &C  ~   <fC      is  the  positive  definite  operator  V       we  have 
that 

<7C°>*  <M>  (2.3) 

for  every  function  ^r  in  the  domain  of  <rC     . 

If  we  order  the  eigenvalues  of  aV       and  aC      lying  below  their 
respective  first  limit  points,  so  as  to  form  increasing  sequences 


and 


£;*£"/*-••,  (2.4) 


E,    *   ££    *     '  '  '  ,  (2.5) 


2 
The  major  part  of  this  Section  and  Sections  III,  IV,  and  V  are 

condensations  of  the  original  development  given  in  [6]  and  [7]  . 


(2.3)  is  sufficient  to  guarantee  [8  ] the  inequalities 

1     ~    ^/   >    C^ '•  /,*,•  •  •  )  .     (2.6) 

We  now  construct  a  sequence  of  Intermediate  HamfltonFans 
satisfying  the  inequalities 

<x'>«<M">±<&—>±<ge>t    (2.7) 

so  that  their  ordered  eigenvalues  satisfy 

E. °i±£j    ±£~-  *±   £i     ,Cf-, *,...)     (2.8) 

and  give  lower  bounds  improving  toward  the  eigenvalues  of  fr  To 

construct  the  ^^we  introduce  a  projection  operator   ^"defined 
in  a  vector  space   S  characterized  by  the  metric  operator  |/ .   The 
inner  product  in   5   is  therefore  g ,  ven  by  <  J  ////#>  .   However> 
elements  of   S  are  taken  to  be  normalized  in  the  sense 

<  J  I  <&)     =  /   .  (2.9) 

In  Appendix  I  we  show  that  the  projection  operator  for  the  sub- 
space  spanned  by  the  elements  ?fe  ,,/,*,--.,  ^  of  a  set  fa] 
linearly  independent  in  S      is 


(2.10) 


where  A       is  the  inverse  of  the  nth  order  matrix  with  element: 
(fc'''lflj)'         It  is  readily  verified  that 


JMatri 
as  a  super 

ad  j  i 


>rscriDtrtoindiCated   bY   b°,d   Synb°,S-      T^   symbol      /applied 
-rscr.pt    to   an   operator   or   matrix   denotes    its    Hermi t ian   adjoint . 


(OV*  =  o~  , 

(2.H) 

0-V  =  omtv  om    =  ^~ 


so  that  the  operator  /^  ^  is  self-adjoint. 

Further,  from  (2.3)  and  Bessel's  inequality  in   0  , 

0      ±     (I/O")    ±     <l/0""}±    <V)     .  (2.12) 

The  relation  (2.12)  forms  the  basis  for  construction  of  the  /^ 
for,  if  we  define  <TL       by 


»m=  76°  -r  I/O",  (2.i3) 


the  inequalities  (2.7)  are  insured  by  (2.12). 

Thus,  provided  the  eigenvalue  equations  of  the  operators  (2.13) 
can  be  solved,  we  have  a  means  of  obtaining  a  sequence  of  lower  bounds 
which  increases  toward  the  eigenvalues  of  <rC 

A  formal  solution  has  been  given  by  Bazley  and  Fox  [6]  which  we 
repeat  in  part  here.   The  eigenvalue  equation  for  the  (TC      , 

may  be  written  with  the  aid  of  (2 . 1 0)  as 

(7i'-E")4m  =  -vom4m 

where 


If  E  is  not  an  eigenvalue  to  ?£  we  may  multiply  (2.15)  by 
the  resolvent  operator  (ft0-  Ej  to  obtain  an  expression  for  the 
eigenf unct ion 


4m  =  -  (H'-Em)"v  EC, 


(2.17) 


Unfortunately,  the  resolvent  operator  is  seldom  known  in  closed 
form  being  generally  only  expressable  as 

(3£*-Em)~'  =  E  /  ^zx+zi 

^*  £■  (2.18)   • 

J*:         £   -  Em 

where  the  4>±    represent  the  discrete  eigenf unct ions  of  'H  °  and 
the  <P  (£) continuum  eigenf unctions  normalized  on  the  energy  scale. 

^   denotes  the  first  limit  point  of  74     ,  which  is,  incidentally, 
also  the  first  limit  point  of  the  PC"*   [6] 

(fC       and  <f(    may  have  common  eigenvalues  in  which  case  the 
solutions  take  a  different  form  [6]  but  the  difficulties  associated 
with  the  resolvent  remain. 

In  the  following  sections  we  discuss  procedures  which  remove 
these  difficulties  and  lead  to  practical  methods  for  computing  lower 
bounds.   Prior  to  this,  however,  we  wish  to  comment  on  the  convergence 
of  the  eigenvalues  of  the  #""  to  those  of  Jr^   . 

It  has  been  shown  that  [6]  ,  provided  H°  and  #  have  completely 
continuous  inverses,  V    is  bounded  relative  to  #*,  and  the  set  {fa} 


is  complete  in   O  ,  the  eigenvalues  and  eigenf unctions  of  the  #t, 
converge  to  those  of  ^C     .   Of  equal  interest  in  applications  is 
the  rapidity  of  convergence.   In  this  connection,  we  demonstrate 
that  the  error  in  a  lower  bound  £.       is  of  second  order  in  the  error 
in  the  eigenf unction  H*      and  in  the  error  in  the  operator  <^£ 
To  this  end  we  introduce  the  orthogonal  complement  t         to  &     de- 
fined by 

l/P"=     ]/'     VOmt  (2.19) 

having  the  properties 

(P"f-  Pm, 
P"jy    =   p~'y  p"  =  yp~  (2.20) 

OmtVpM=  PmtVCT=  0  . 

We  may   now  write    J~£    as 

H  -   H°  t  y( 0"  +  P") 

(2.21) 

=   Hm   +  V?"   > 

and  one  of  its  eigenf unctions  as  the  sum 

4  -  4m  ■>  <f,  (2-22) 

where 

(4  /«(>     =    1  ,         <JlS>    --0  .  (2.23) 

The  energy     £     corresponding    to     dp     is    then 

l  -n**    ( * 

We  point  out    that    /        <fJ      does    not    in  general    vanish. 


(2.24) 

=  £""<V7  «0  *  <<sig€-is>  +  <4ivpml4>. 

If  we  denote  r  ^i  by  O   and  retain  only  terms  of  second  order 
we  find  for  the  error  in  the  lower  bound 

04  A~=    £  -£" 

(2.25) 

=  <s/%~/s>  -E~au>  t  <j'/i//j'>  , 

observing  that  the  first  two  terms  arise  from  the  error  in   <f*   while 
the  last  arises  from  the  error  in  (TC 

We  give  for  comparison  the  well-known  equivalent  relation  fcr  the 
error  in  a  upper  bound  obtained  by  the  variational  principle.   V/e 
take   5P  to  be  a  normalized  trial  function  orthogonal  to  all  eigen- 
functions  of  </v  with  energies  less  than  E    .      We  wr i te   J  as 

<£>    =    4$    +    Jf  (2.26) 

where 

<$/$>=  i       ,    <4I%>  =  o  .        <2-"> 

If  W     is  the  upper  bound  obtained  with  ^  .i.e., 

A/  -    <<£IMI  $> 

(2.28) 

=    E<4IJ>   +  <*/X/f>     , 

we  find  for  the  error  in  the  upper  bound 


Upon  comparing  (2.25)  and  (2.29)  we  see  that  the  first  two  terms 
of  (2.25)  are  of  the  same  general  form  as  the  terms  of  (2.29)  and   if 

<£*   and   y  do  not  differ  greatly  it  is  reasonable  to  suppose 
that  these  terms  will  be  comparable  in  magnitude.   If  this  is  the 
case  we  see  that  the  error  in  the  lower  bound  will  be  greater  by 
approximately 

A'- A*   e*    <<T7///<T>.  (2.30) 

We  may  therefore  expect  that  procedures  based  on  the  intermediate 
Hamiltonians  given  here  will  not  demonstrate  convergence  properties 
comparable  to  upper  bound  procedures  based  on  the  variation  prin- 
ciple.  We  find  confirmation  in  our  application  to  the  helium  atom 
in  Section  IX. 


III.   THE  METHOD  OF  SPECIAL  CHOICE 

Bazley  [5]  proposed  the  special  choice  of  the   y*A 

/•*   ~  P"'    *&*         ,       ( *    :/,*>■■"+)>       (3-0 

where  the  5^  are  again  elgenf unctions  belonging  to  the  discrete 
spectrum  of  rC     .   With  this  choice  (2.17)  becomes 


-  -Z4z 


(3.2) 

C+  +2 


Upon  substitution  of  (3.2)  into  the  eigenvalue  equation  for  <rC      , 


we  f  i  nd 


or  (3.3) 

EC,  ¥>'  =-o"^E c9%  . 

By  forming  the  Inner  product  of  (3.3)  successively  wi  th  a*  a.*,  •  •  ',  *,* 
we  obtain  the  /n   equations 

10 


11 


which  constitute  a  matrix  eigenvalue  equation  of  order  fr\       .   This 
equation  can  be  put  in  more  conventional  form  if  we  successively 
multiply  on  the  left  by  the  diagonal  matrix  with  elements  ( £ *'  ~ £■    '"*} 
and  the  inverse  of  the  matrix  with  elements 

<fji/if1y  =  <4>°jv"i<p;>  (3-5) 

which  we  have  previously  designated  by  f\      .      The  resulting  eigen- 
value equation 

£[(E}-Em)S,i+A"J=  o  ,         (3-« 

differs    from  a  Ray le igh-R i tz   calculation   based   on    the    trial    function 

$=IlAj+l  (3-7) 

4.  =•/ 

only  in  that  the  matrix  |r   is  replaced  by  f\ 

(3.6)  determines  /yi      of  the  eigenvalues  of  *%.      and  the  cor- 
responding eigenf unctions  are  of  the  form 

is  i 
The  remaining  eigenf  unctions  of  rt      are  all  eigenf  unctions  of  *[, 
different  from  $,,  ¥4   ,*''     ¥"*    •   These  have  corresponding 
eigenvalues ,  i.e., 

%"  +  •  =£;*;  ,(/>„).     (3.9) 

When  these  pers  i  s  tent  eigenvalues  are  ordered  with  those  determined 
from  (3.6)  in  a  non-decreasing  sequence 


12 


E~*E"±  •  •  .  ,  (3.10) 

we  have  the  lower  bounds 

e;  -  e1   ,  </-<*,...)  .    (3.10 

This  procedure  is  sometimes  simpler  to  apply  from  a  computational 
viewpoint  than  the  other  procedures  covered  in  this  paper  and  can  be 
expected  to  give  relatively  good  estimates  for  excited  states  much 
more  easily.   We  note,  however,  that  for  problems  where  a  portion  of 
the  spectrum  of  &      is  continuous,  the  ^  given  by  (3.1)  do  not 
constitute  a  complete  set,  and  the  procedure  cannot  converge. 

For  an  application  to  a  problem  involving  continuous  spectra 
we  refer  to  Bazley's  calculation  on  helium  [5  ].   For  an  application 
to  a  problem  where  ctt      has  only  bound  states  see  Section  VIM. 

This  procedure  has  been  generalized  [6,7]  to  include  ^  given 


by 


1*       ^   ft     fy    +)' 


(3. 12) 


The  eigenvalue  equation  equivalent  to  (3.6)  is  slightly  more  in- 
volved but  may  lead  to  simpler  calculations  in  some  cases  [7  ]  . 


operator 


IV.   THE  METHOD  OF  TRUNCATIONS 

This  method    [6,7]    involves    replacing   JC      by    the   smaller 
tor5 

(4.1) 

+  E°    [  l  -  El+lK+ll] 

to  obtain  the  truncated  intermediate  Hami 1 tonians 

%'•"  =  ?{*''+  vo"  .  (!t.2) 

These  satisfy  the  inequalities 

<7t£'">±<x'"''">4<M">±<ft>    (4.3) 

so  that  their  eigenvalues  give  lower  bounds 

E*f  *  £/,   </?<*,-).      <-» 

The  advantage  of  this  procedure  is  that  the  resolvent 
\  oi,         ~~ C-  /      is  known  in  closed  form, 

/  -  27/00;/  d.5) 


+ 


E1+*  -E 


£,  sr< 


The  £.  »    are   taken  to  be  ordered  according  to  (2.k) . 

13 


14 


The  relation  equivalent  to  (2.17)  is 
£ 


which,  when  substituted  into  the  eigenvalue  equation  for  $-£    ', 
leads  to 


Ct.6) 


r 


.'Jr   <<t>llvl*<>    jl. 


e:„  -  e 


£,  *n 


4>:  + 


C,   =o. 


(4.7) 


By  taking  the  inner  product  of  (4.7)  with  the  elements  stj*      wh.  * 

(^  *    J    7^  Jt  i 
*        .      th 
'''if*   we  obtain  an  n   order  matrix  eigenvalue  equation 


E 


El  -  E'-~ 


<Vfitivt>,.y  -  Z^ivi+tX+llvitJ 


e;.,  -  e 


£,/n 


C<  =  o, 


(4.8) 


(f  =  '**,'-'  ,m)     > 


whose  solutions  give  all  elgenf unctions  and  eigenvalues  of  the  ?€ 


C 


different  from  those  of  /t   .   If  there  are  less  than  I  +  /»   such 

solutions,  it  remains  to  be  determined  whether  any  of  the  eigenvalues 

E  ,  ,   £ i   ,'"t   t- /»  persist  in  the  spectrum  of  f"C     '        .   Procedures 

for  finding  such  eigenvalues  and  their  eigenf unctions  are  given  [6,7] 

but  we  do  not  cover  them  here  as  they  are  not  needed  in  our  application 

of  this  method.   Finally,  it  is  shown  [6,7]  that  all  other  eigen- 

functions  of  a~C        belong  to  the  infinitely  degenerate  persistent 
e  • 

eigenvalue  *-  g  +  t    . 

Equation  (2.25)  for  the  error  in  a  lower  bound  requires  modifi- 
cation for  this  procedure.   We  have,  instead  of  (2.21) 

so  that  Zj       becomes, 

o*  A'  =  £  -£*"•    - 

=  <S I &'"!£>  -El'"'<SlSyf<S'IVIS'>  (*-,0) 

The  additional  term  cannot,  in  general,  become  arbitrarily  small  unless 
vL         has  bound  states  only,  showing  that  the  procedure  is  not  con- 
vergent for  problems  involving  continuous  spectra. 

Nevertheless  we  note  that  this  procedure  can  always  give  results 
superior  to  those  obtained  by  the  method  of  special  choice  for  it 
contains  the  method  as  a  special  case.   This  is  seen  by  observing 
that  when  the  AiA    are  chosen  by  (3.1),  (^.8)  reduces  to  (3-^). 

Bazley  and  Fox  have  applied  this  procedure  to  helium.  We  give 
some  additional  calculations  in  Section  IX. 


V.   A  SECOND  METHOD  OF  TRUNCATIONS 

Bazley  and  Fox  have  given  another  method  of  truncation  [6,7] 
which  can  be  convergent  for  problems  involving  continuous  spectra. 
To  illustrate  this  method  we  write  <rC     as 

X   =    7t£'°  +(X-  2t*'°)  ■  (5.0 

where  the  <rC         are  given  by  (4.1).   Since  the  operator  '  "I "  *"        ' 


is  positive  the  operators 


(5.2) 


are  intermediate  between  <r£    and  <r£  .    is        is  given  by  (2.10) 

with  V      replaced  by  (  <K  ~  at      J. 

1/  ^j/* 

The  solutions  to  the  eigenvalue  equations  for  the  vl>  are 


developed  exactly  as  in  Section  IV.   The  equations  corresponding  to 


(k.S)    and    (4.8)    are 


4Mtm-E 


<  +:  1 1//&> 


p   <  <r*i  rip*/     <£"     + 
-      £1    -  £*" 


(5.3) 


(n-H9>*)fi* -£.<+'*  !"*>+* 


L  ^  f  /  E- 


£j/r* 


C, 


and 


To  distinguish  the  two  truncation  procedures  we  will  hereafter 
denote  the  present  procedure  by  truncations  ii  and  the  previous  pro- 
cedure by  truncations  J. 


16 
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z 


<r,ix-H£t'i^>  + 


(5.4) 


<(X-X")Ml(X-X**)*>-£<A>M+2<#W*t 

-7^7- 

( /  -  ***»'  "  j  "*'  . 

Procedures  for  finding  persistent  eigenvalues  parallel  those 

r° 

for  the  method  of  truncations  i.   Here  E-£+i    is  also  an  eigen- 
value of  infinite  multiplicity. 

Sufficient  conditions  for  convergence  for  this  procedure  are 
just  those  stated  for  the  vC      in  Section  II  with  r        replaced 
by  \H~7i      J,    and  f{,      replaced  by  f\,     '  We  point  out  that, 

should  tfT   have  a  spectrum  extending  to  infinity,  then  (&C~?{       / 
is  not  bounded  relative  to  <fl         for  £      finite. 

This  procedure  has  been  applied  to  the  helium  atom.   Our 
results  are  given  in  Section  IX.  We  find  the  lower  bounds  given 
by  this  procedure  to  be  inferior  to  those  of  the  other  three 
procedures  covered  in  this  paper.   For  a  possible  explanation 
we  refer  to  Section  VII. 


VI.   ELEMENTS  CONTAINING  THE  INVERSE  OF  THE  RESOLVENT 

We  present  here  a  new  procedure  for  solving  the  eigenvalue 
equation  of  the  intermediate  Hamiltonians  (2.13).   Our  procedure 
avoids  the  appearance  of  the  resolvent  operator  (2.18)  by  selecting 
the  *ft.i,     according  to 

r,  =    V(9tm-Em)fi  ,  6=/*,  ••  ),».» 

where  the  ti     are  essentially  an  arbitrary  set  of  elements  in  the 
doma  i  n  of  <rC    . 

In  the  development  of  this  procedure  it  is  convenient  to  introduce 
a  family  of  ^A^    more  general  than  these  of  (6.1)  defined  by 

fut*)   =  V"(k'-e.)h  ,  (**i,z,-),  <6-2> 

where  €      is  regarded  initially  as  an  arbitrary  scalar  parameter. 
We  require  of  the  T^     that  they  be  linearly  independent  and  such 
that  the  jh±  (£)  are  square- i  ntegrable  in   ^   ,  i.e.,  such  that 
\  ft*  \€  )lrlfii(€)/         isfinite.   We  do  not  requ  ire 
that  the  operator  r     (r(~€/   considered  as  a  linear  transfor- 
mation from  the  set  llil    to  the  set  ifli)    be  invariably  non- 
singular  for  all  values  of  £ 

The  author  Is  indebted  to  his  advisor,  Professor  P.O.  Ltiwdin 
for  suggesting  this  choice. 
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From    (2.10)   we  obtain  for    the  operators   O (€  )    projecting  on 
the  subspace  spanned   by     *P  1  (€)  }   4i  &   (£)  j    '    '    '    ,    Jo,*,  (£  )  t 

0(€)    =  IZ1  //*+<€)>  Ay <^/C€)/l/.  (6.3) 

Using  these  we  construct,  as  in  (2.13),  a  family  of  intermediate 
Hami 1 tonians  parameterized  by  fe  > 

7€"ce)  =  9£*+  l/0"(6)  .  (6.4) 

According    to    (2.17)    their   eigenvalue  equations, 

H  (C)  <p  CO   =   E    <PC€)  ,  (6.5) 

are  satisfied,  for  eigenvalues  different  from  the  £_    ■   ,  by  eigen- 
functions  of  the  form 

4%)  =  (9£'-Em)"l/EC+fr<6),     (6.6) 

with  the   C*  given  by 

Cx    =  £!  A+j(€)</if(ejl VI <£(£))  .     (6.7) 

Using  (6.2)  in  (6.6)  gives  an  expression  for  the  eigenf unctions  in 
terms  of  the  f-      , 

J7e)  =  (K*-E~)~'£&(X'-e)&  .     «•«) 


We  inquire  now  as  to  whether  there  exist  specific  operators 

/»V  *•«•! 


<rl    ' £   /of  the  family  (6.^)  having  the  eigenvalue  E-      for,  in  th 
event,  the  eigenf unction  (6.8)  reduces  to  the  finite  linear  combi- 
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nation  of  the 


h 


4?e")  =  £C*L 


(6.9) 


If  such  operators  exist  they  will  satisfy  eigenvalue  equations  of 
the  form 


(W°-  E)£Cj.h   =  VO«\£~i  Z  C+  £ 


/n 


-0"f(E")Vi:C.fx 


(6.10) 


*,r~ 


By  forming  the  inner  product  of  (6.10)  successively  with  ^ft,  ( £.       ) 
,4& f    ( E.      Jj  '   '    'j  Ax^KC.    /we  obtain  the  matrix  eigenvalue  equation 


(6.11) 


E [  <(K'-E~)fjlV-'l(M'-Er)PA> 
<(9e'-E~)PilV"l(VC-E~)fi.> 


which  has  solutions  provided 

r 


Jet 


+  <hl9f-E'"lh> 


(6.12) 


=  0 


for  real  £L       .   The  existence  of  real  roots  to  the  secular  equation 
(6.12)  is  thus  the  criterion  for  the  existence  of  operators  at     \  L.     / 


.    We  observe  that,  since  the  resolvent  has  been  eliminated  from 
in  going  from  (6.8)  to  (6.9),  (6.9)  is  not  restricted  to 
eigenvalues  different  from  the  £*      . 


4' 
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The  left  side  of  (6.12)  Is  a  polynomial  of  degree  Z/r\  in  £.    and  has 
an  even  number  /Vh  of  real  roots  (O    4z /rnl£:  £ ^ )  .      Each  such  root  £" 
is  an  eigenvalue  of  an  intermediate  Hamiltonian  oC     (Ej   /•   If 
there  are  M.       distinct  roots  (4.  — ""J,   we  find  from  (6.12)  eigen- 
values to  the  A     operators  rC      ( £.,    J    ,  #T  ( £ '      )       '    '    '  f 
a*"      **-j£/      Should  a  root  £.4    be  multiple  it  is  a  degenerate 
eigenvalue  of  aC   (c;/   The  associated  eigenfunctions  <%J     (p'") 
<t>   (£*)/  '  *j  *T  (£7)  are   determined  from  (6.11). 

The  fact  that  the  procedure  supplies  only  one  eigenvalue  for 
a  given  operator  is  a  serious  disadvantage  for  without  the  remaining 
eigenvalues  we  cannot,  in  general,  locate  the  known  eigenvalue  £_  ^ 
in  the  spectrum  of  o(    (  k 4 ) and  are  consequently  unable  to  deter- 
mine to  which  eigenvalue  of  ^"£   it  is  a  lower  bound.   However, 
from  the  relation 

<X'>*<  2€7£p>,  (6.i3) 

9 

£1  *  £7  ,  (* ='.*,  ■■)  ,     <6-'« 

which   serves    to  prove    the  following   Theorem. 

THEOREM:       If    the   known   eigenvalue    £  ;of  of(Ej  ^satisfies 

E'm  *  £?*  £«.,     ,  (6.-5) 

The  subscript  on  £-jL      in  (6.14)  enumerates  the  ordered  eigen- 
values of  2s"}rm)   • 


we  have 
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and  if  upper  bounds  to  oc  -  /  eigenvalues  of  {^   lie  below  C-  j  , 
then  £  J   is  a  lower  bound  to  £.oc  .   For  the  ground  state  the  theorem 
reduces  simply  to 

E^    *£l  .  (6-16) 

Fortunately,  (6.16)  is  sufficient  to  identify  lower  bounds  to  the 
ground  states  of  many  atomic  and  molecular  systems. 

In  connection  with  the  convergence  of  this  procedure  we  first 
consider  the  entire  family  of  operators  (6.4)  without  regard  to 
whether  they  have  solutions  determined  by  (6.11).   Since  these  are 
just  special  cases  of  the  intermediate  Hamiltonians  of  Section  II, 
the  conditions  stated  for  convergence  there  are   sufficient  for 
convergence  of  any  member  of  (6.4)  to  yv      .   However,  the  com- 
pleteness of  the  elements  (6.2)  requires,  in  addition  to  the  com- 
pleteness of  the  Fx     ,  that  the  transformation  defined  by  (6.2) 
be  non-singular.   In  Appendix  I  we  show  that  the  necessary  and 
sufficient  condition  for  the  elements  At/  (6),  4±z(£  h  '  '  ')fi#>'£' 
to  be  linearly  dependent  is  that  some  ei  genf  unction   tJ     of  /&■ 
be  i  nc  1  uded  i  n  the  span  of  fi  ,fg  ,*  •  •  t  flu     and  tha  t     ^    ^  £j    ' 
Therefore,  if  we  separate  the  family  (6.4)  into  those  operators 
for  which  €   7^  £ X  j   4s'***'  ••  ,    and  those  for  which 

£    =    C.  ^     ,  X'  *  ft^t'*    '    ,  the  former  will  converge  to  rC     as 
the  TjC     become  complete  while  the  latter  will  not  since  fC    (£j'   s 
will  agree  with  Tv      on    '4     >  '«e»i 

WCE})*}  =  £}  +;.  (6. t7) 

This  fol  lows  from  the  eas  i  ly  deduced  resul  t  that  if  / ;  ,fx  ,  •  -  -;  f/m 
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contains  $: 

<+;/v07£;)i+p  =o  .  (6. is) 

We  return  now  to  the  members  of  (6.^)  which  have  an  eigenvalue 
and  eigenf unction  determined  by  (6.11).   We  emphasize  two  aspects 
of  the  solutions  to  (6.11)  in  which  they  differ  from  the  solutions 
for  the  other  procedures  covered  in  this  paper.   First,  it  is 
entirely  possible  that  a  given  equation  will  have  no  solutions  at 
all  corresponding  to  a  situation  in  which  (6.12)  has  no  real  roots. 
Second,  we  have  no  guarantee  that  the  subspace  of  is         l  t.         J 
will  contain  that  of  Cs     (£   /s  i  nee,  in  general  ,  A  v  ^  ci   v  . 

r-  •»"  r-  ■*»  +  ( 

We  therefore  have  no  assurance  that  the  sequence  £■  V ;  c.    *  '   .   *  *  " 
is  monotonic  increasing.   Because  of  these  complicating  features 
we  discuss  only  a  special  case  which,  however,  is  the  one  of  in- 
terest in  applications.  We  suppose  that,  by  a  judicious  choice 
of  the  //.       ,  we  have  succeeded  in  obtaining  a  set  of  equations 
(6.11)  of  orders  Q,    Q+J   ,    4*2   ,       etc-  whose  secular  equations 
have  at  least  one  pair  of  real  roots.   Further  we  suppose  that 
each  member  of  the  sequence  /_  >  ,  £-    *     ,    '    '   '      formed  of  the 
•4 ty     roots  of  the  secular  equations  satisfies  the  theorem  pre- 
viously proved  for  some  value  of  oc      .  We  then  have 

We  now  make  use  of  the  convergence  properties  established  above 
to  state  that  as  /rj       increases  without  limit  and  the  r*c 
become  complete  one  or  the  other  of  the  equalities  in  (6.19)  will 


2k 


hold  so  that  sequence  will  (on  the  average)  either  decrease  to- 
ward u^  oc   or  increase  toward  £-<k.      We  are,  of  course,  in- 
terested only  in  lower  bounds  belonging  to  sequences  exhibiting 
the  latter  behavior. 

As  with  the  method  of  truncations  i,  this  procedure  contains 
the  method  of  special  choice  as  a  special  case.   When  the  /^  are 
chosen  as  eigenf  unctions    '/*- of  &£     ,  (6.11)  may  be  seen  to  reduce 
to  (3.*0. 

We  mention  a  feature  of  this  procedure  which  is,  to  our  knowl- 
edge,  unique  among  procedures  based  on  intermediate  problems.   It 
is  that  the  procedure,  while  requiring  a  knowledge  of  the  eigen- 
values of  y^C      »  does  not  use  its  eigenf unctions .   This  can  be  a 

definite  advantage  in  applications  to  systems  such  as  diatomic 

r~  ° 
molecules  where  the  £_  *     mav  be  accurately  determined  but  where 

rh>°  10 

the    i   j    are  n°t  known  in  closed  form. 

We  are  of  the  opinion  that  the  success  of  this  procedure  in 

our  application  to  helium  in  Section  IX  lies  in  the  ease  with  which 

■^kj.   may  be  chosen  which  lead  to  good  intermediate  eigenf  unctions . 

This  feature  is  shared  in  large  measure  by  the  method  of  truncations 

?.   We  close  the  theoretical  part  of  this  paper  with  some  comments 

in  this  connection. 


10  C  °  eb° 

For  diatomic  molecules  the  t-^   and  i/,   are  constructed  from 

the  solutions  to  the  hydrogen-molecule  ion  Hamiltonlan  whose  eigen- 


functions  are  not  known  in  closed  f orm  [  9  ] 


VII.   DESIRABLE  FORMS  FOR  INTERMEDIATE  E IGENFUNCT IONS 

We  consider  the  error  in  a  lower  bound  due  to  the  error  in 
the  eigenf unction  which  according  to  (2.25)  is  \ol  1H>      ~ £.       I O/ 
In  general  we  expect  a  reduction  in  O       to  improve  the  lower  bound 
so  we  wish  to  have  the  quanti  ty  I  \  ty      I  *r  '     I     large.   This  is 
most  easily  accomplished  by  selecting  the  A) ^     so  that  their  con- 
tribution to  the  eigenf unction  is  a  linear  combination  of  functions 
known  to  serve  as  a  good  trial  function  for  a  variational  upper 
bound  calculation.   If,  e.g., 


$*E/I;fi  (7.0 


4=1 


is  such  a  good  trial  function  we  would  choose  jb^     given  by 

^x  =     Y~' Ta         ,  (  j.=    *,*>-,">   )   (7.2) 
for  the  method  of  truncations  i  ,  and  by 

fi      =  V"(%'-£~)f,  ■,(*' l.t.  ■■;*)     (7.3) 

for  the  method  of  Section  VI.   The  respective  eigenf unctions  are 
then  of  the  desired  form, 


This  is  invariably  true  for  the  ground  state. 
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-f  I 


**'"*  -£[£(£: -l?*i(E;,r£*~) <r'lh>  *: 


and 


+ 


A 


c. 


(7.4) 


^?£*j  -£&/, 


(7.5) 


However,  if  we  attempt  a  similar  choice  for  the  method  of 
truncations  i  i  ,  i.e., 

ft     =  (9t-9i''')"fi   ,(j.t.*,...,m),    (7-6) 

we  find  that  we  are  required,  in  evaluating  the  elements  of  (5.*0, 
to  compute  the  inner  products, 


</>*  m-  n  '••/?,>  =  <fii(7t-9t''*)itj>.  (7.7) 

This  is  generally  not  practical,  being  a  task  of  the  same 
nature  and  difficulty  as  is  found  in  applying  the  resolvent  operator 
in  Section  II.  We  are  of  the  opinion  that  our  poor  results  for 
helium  given  in  Section  IX  are  a  consequence  of  our  not  being  able 
to  select  our  -4lj  accord i ng  to  (7-6). 


VIII.   A  TWO-PART  I  CLE'  HARMON  I C  WELL 

We  determine  lower  bounds  to  the  eigenvalues  of  a  system 

consisting  of  two  electrons  confined  by  an  isotropic  harmonic 

12 
well,  having  the  Hamiltonian 

7c-iv.*-ir:*-k "V^~v^ .  (8.,» 

Def  i  n  i  ng 

and 

1/  =     —  (8.3) 

we  wr  i  te  v~t     as 

X  =  M"  +  V  ,  (8.4) 

and  see  that  it  is  of  a  suitable  form  since  JV      and  /r   are  positive 
definite,  and  the  ei genf unctions  and  eigenvalues  of  ^C     may  be  de- 
termined from  the  known  ei genf unctions  and  eigenvalues  of  the  one- 
particle  Hamiltonian 

X,  -  -kF*  +  i ■«**/".  (8.5) 


1  2 

Atomic  units  are  used  throughout  Sections  VIII  and  IX, 
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We  use  the  method  of  choice  of  elements  detailed  in  Section  III. 
We  note  that  the  Hamiltonian  (8.1)  is  well  suited  to  this  method  since 
its  spectrum  consists  entirely  of  bound  levels. 

We  outline  here  the  procedure  for  constructing  the  ei genf unctions 
and  eigenvalues  of  At  from  those  of  tfT,  .   Since  SV  °  may  be 
written  as  the  sum 

7i      =  7i,(i)   +  ?4,(2)  ,  (8.6) 

0 

the  eigenfunctions  of  SV   are  linear  combinations  of  products  of 
the  eigenfunctions  of  Yl,(/)   and  7i,  (2) ,  and  its  eigenvalues  are 
sums  of  their  eigenvalues.   It  is  shown  [10,  p.  1663]  that  the 
eigenfunctions  of  #V,  are 


*  ,  (8-7) 

i—i*-i 

and  that  the  corresponding  eigenvalues  form  a  discrete  spectrum 
extending  to  infinity, 

r,,m    =f»,+±)cu*,      (*,=  /,*,-  ■■    ),        (8-8) 
with  all  levels  above  the  lowest  degenerate.   In  (8.7)  the  )//»nare 

/  *  +  i 

spherical  harmonics  and  the  [_,  ,  ,  f       /   are  Laguerre  polynomials 

tm  "J* '  1 

[10,  p.  78^+]  .  Jf      and/tM  are  the  usual  total  and  magnetic  orbital 

angular  momentum  quantum  numbers  and  /r,    is  a  radial  quantum  number. 
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In  our  calculations  we  consider  only  singlet  S  states  of  crX 
so  that  we  need  only  singlet  S  states  of  aC  .  The  space  parts 
of  these  are  linear  combinations  of  products  of  the  form 

^/,/m,>*  CD    y.£t/n',/»' '(*)    ,  f/m'--/n»)    (8'9) 

symmetric  in  the  coordinates.   For  convenience  we  designate  them 

2 
with  the  usual  spectroscopic  notation,   (is  ),  (2p  *+p)  ,   etc. 

Hylleraas  [11]  has  shown  that  these  space  parts  may  be  written  as 


Or,  «  „  a)    =  ^   [  S^jCr.)    S<r>,t(rj)]  Y^oCff^), 


(8.10) 


(8.11) 


when  /r»'^  sr\    .        O/rt   £    's  tne  normalized  radial  part  of  (8.7)  and 

(s/%     is  the  angle  between^"  and  /^  .   The  symbol  Q      here  represents 
one  of  s ,  p ,  d ,  etc . . 

In  this  notation  the  eigenvalue  equation  for  wC 


o 
is 


X'Csrifsn'a)      =     (/»  +syt'+  /)co2(/»Qsr,'q)    .    (8.12) 

We  list  for  reference  the  seven  lowest  singlet   S   eigenf unctions 
of  fC       and  the  corresponding  energies: 


State  Pes  ignation  Energy 

4>°,                                                       (Is2)  3  LU* 

+\                                                   (ls3s)  5u>2 
+  °3                                                       (2P2) 
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State  Pes  ignation  Energy 

f°+                                                     (ls5s)  jcu2 

4>Z  (3s2) 

+  1                                                      (2p4p)  // 

+  °r                                                      (3d2)  // 


Having   determined    the    lower   part   of    the   spectrum  of  <fu      we 
define    the  sfl    ■    according    to    (3.'),    i.e. 

pi     -  1*    +1        ,      (*'  *   'j*,-"*T),       (8.13) 

and  compute  the  elements 

The  procedure  used  to  evaluate  the  integrals  is  covered  in  Appendix  I 

For  an  n   order  ca  1  cu  1  at  ion  ( *r\£   "we  select  a  submatrix  of 
order  /*J   from  the  matrix  formed  by  the  elements  (8.12)  and  invert 
it  to  obtain  the  matrix  /%      .   Having  determined  the  matrix   XI 
we  form  the  eigenvalue  equation  (3.6) 

E[(E;  -E~)S*j  t  A^JC^o,      (8.15) 

which  may  be  solved  for  the  eigenvalues  £      and  the  eigenvectors 
W    by  conventional  methods.   When  the  eigenvalues  of  (8.15)  are 
ordered  with  the  eigenvalues  which  persist  from  <rC     we  have  lower 
bounds  to  the  eigenvalues  of  qC 

In  order  to  get  an  idea  of  the  accuracy  of  the  lower  bounds  we 
have  computed  upper  bounds  by  the  well-known  Ray  lei gh-R i tz  procedure 
based  on  trial  functions  of  the  form 
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$ =  c f\x  +1    c„±  7) . 


(8.16) 


The  upper  bounds  are  given  as  the  solutions  of  a  matrix  equation 
identical  to  (8.13)  except  that  f\       is  replaced  by  V       with 
elements  \T+lpl   r*/  ,    i.e., 


EKE',  -E")J,u  *  VfAAx  =  o  , 


<V  =  '*Z, 


For  our  calculations  we  chose  the  value  0.64  for  cu 
choice  gives  to  the  ratio 


ff  = 


<4.ivi4.> 


<4.iffl4.> 


(8.17) 


This 


(8.18) 


a  value  of  about  0.3  which  is  comparable  to  that  for  the  helium  atom. 

Our  results  are  given  in  Tables  1,  2,  and  3.   Table  1  displays 
the  matrices  V  V      /.its  inverse  f\       ,    and  the  matrix  \f       used  in 
the  upper  bound  calculations.   It  is  interesting  to  note  the  close 
correspondence  between  /\      and  \f        .   Our  upper  and  lower  bounds 
are  compared  in  Table  2.   The  first  order  calculations  were  based 
on   r  /  ,  and  the  third  order  on   ^/  ,  i  *    and  13  .      The 
elements  added  for  each  of  the  fourth  order  calculations  are  given 
in  the  table.  We  compare  the  normalized  eigenf unctions  obtained 
from  the  seventh  order  upper  and  lower  bound  calculations  in 
Table  3.   We  find  that  eigenf unctions  representing  the  same  states 
are  substantially  the  same.   Outstanding  is  the  nearly  complete 
agreement  of  the  eigenf  unctions  for  C.  q.    • 

We  see  that  the  rate  of  convergence  of  the  upper  and  lower 
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Table  2.   Upper  and  lower  bounds  to  the  singlet  S  states  of  the  harmonic 
wel  1  . 


Order 

El 

E2       h 

\ 

h 

h 

E7 

Upper  Bounds 

1st 

2.5583 

3rd 

2.5030 

3 

7872   3.8383 

4th  (ls5s)  2.5023  3.7729  3.8191  4.9882 

4th  (3s2)  2.5019  3.7713  3.8370  4.9339 

4th  (2p4p)  2.4998  3-7644  3.8374  5.0508 

4th  (3d2)  2.5021  3.7670  3- 8159  5.0428 


7th         2.4968   3.7290   3.7824   4.8201  5.0100  5.0673  5.1156 


Lower  Bounds 
7th         2.4730   3-6430   3.7407   4.7933  4.8010  4.8765  4.9879 


4th  (ls5s)  2.4646  3.6057  3.7152  4.4800  4.4800  4.4800  4.8455 

4th  (3s2)  2.4640  3.6007  3.7129     "      "      "  4.8314 

4th  (2p4p)  2.4684  3.6190  3.6997     "      "      "  4.8646 

4th  (3d2)  2.4652  3.6076  3.7222      "      "      "  4.8667 


3rd         2.4630   3-5995   3.6975     "      "      "    4.4800 
1st         2.4213   3.2000   3.2000      "       "       "       " 
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bounds  are  generally  comparable.   However,  since  we  do  not  know  the 
exact  eigenvalues  we  cannot  make  an  accurate  judgement.   We  find 
nothing  in  these  results  which  contradicts  our  assertion  in  Section 
II  that  the  lower  bounds  give  poorer  estimates  than  the  upper  bounds. 
We  would  not  judge  that  either  the  upper  or  lower  bounds  exhibit 
outstanding  convergence.   We  feel  that  this  may  be  attributed  to 
the  functional  form  of  the  *Ta  .       It  is  known  that  sets  of  functions 
whose  angular  dependence  is  given  by  the  spherical  harmonics  do  not 
exhibit  rapid  convergence  for  operators  containing  the    ~ 
operator  [12]  .   It  is  a  shortcoming  of  the  method  of  choice  of 
elements  that  such  functions  could  not  be  avoided  in  this  calculation. 


IX.   LOWER  BOUNDS  TO  THE  GROUND  STATE  OF  HELIUM 

The  helium  atom  serves  as  an  ideal  system  for  the  testing  of 
lower  bound  procedures.   Its  non-relat i vis  tic  Hamiltonian  possesses 
a  typical  atomic  spectrum  consisting  of  an  infinity  of  bound  states 
lying  below  its  first  limit  point  at  -2  [13]  and  a  continuum  in  which 
are  imbedded  other  bound  states  above  this  point.   At  the  same  time 
the  Hamiltonian  is  of  a  relatively  simple  form  so  that  the  procedures 
do  not  lead  to  unmanagable  calculations.   Further,  as  a  result  of 

highly  accurate  calculations,  the  ground  state  is  known  to  lie  be- 

13 
tween  the  1 imi  ts . 

-2.903726615  (au)  £  £*,  £■  -2.90372^375  (au)  (9.1) 

so  that  the  error  in  a  given  lower  bound  may  be  precisely  determined. 

Bazley  [5]  and  Bazley  and  Fox  [7]  have  previously  applied  the 
method  of  intermediate  problems  to  the  helium  atom.   The  best  estimate 
given  for  the  ground  state  is  -3.0008  (au)  [7]  ,  obtained  as  the 
first  eigenvalue  of  a  truncated  intermediate  Hamiltonian  (rC 
In  this  section  we  give  lower  bounds  to  the  ground  state  using  the 
procedures  of  Sections  IV,  V,  and  VI. 


These  bounds  are  those  given  by  Pekeris  [1^,  15]  •   The  lower 
bound  was  determined  from  Temple's  formula  [16  ] using  the  1078  term 
variational  wave  function  which  gave  the  upper  bound.   A  better  upper 
bound  having  an  estimated  error  of  about  one-half  that  of  Pekeris  has 
been  given  by  Schwartz  [171  •   From  his  work  we  estimate  that  the  error 
in  the  lower  bound  of  (9.0  is  about  1,000  times  greater  than  that  in 
the  upper  bound. 
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We  consider  the  non-relati vis  tic  helium  Hamiltonian, 

U  =  -iV,'-iv/-p-   -F    *  p    >         (9-2) 

jhich  may  be  put  into  the  form  (2.1)  by  taking  7€     as 
md  V       as 

V  ~  —    ■  (9.M 


r. 


M° 


may  be  written  as  the  sum  of  two  independent  hydrogen- 1 i ke 
Hami 1 tonians , 

h  -  -k  v*-y  ,  (95) 

whose  solutions  are  well  known  [18]  .  v\      is  therefore  solvable  for 
its  eigenvalues  are  just  summed  pairs  of  the  eigenvalues  of  fi    and 
its  eigenf unctions  are  products  of  the  corresponding  eigenf unct ions 
properly  symmetrized  to  account  for  the  Paul?  Principle. 

Since  we  estimate  only  the  ground  state  which  is  singlet  S  in 
character,  we  require  only  the  lowest  singlet  S  states  of  ^v  . 
These  consist  of  its  ground  state  and  a  sequence  of  excited  states 
with  eigenvalues  converging  toward  the  first  limit  point  at  -2, 


(9.6) 


39 


The  corresponding  ei genf unctions  have  as  space  parts  the  symmetric 
products 

(9.7) 

where  the  f\j,o  are  the  normalized  radial  wave  functions  of  1*      having 
zero  angular  momentum. 

Having  constructed  the  necessary  solutions  to  £tT  we  turn  to 
a  consideration  of  a  choice  of  the  elements  fl ji    ,   which  are  essentially 
unrestricted  in  the  procedures  we  employ  here.   Our  aim  is  to  select 
the  -itji    so  as  to  lead  to  good  intermediate  wave  functions  as  dis- 
cussed in  Section  VII.   At  the  same  time,  we  wish  to  avoid  elements 
which  lead  to  unduly  complicated  calculations.   Elements  which  meet 
these  requirements  very  well  (for  the  procedures  of  Sections  IV  and 
VI)  are  those  selected  so  as  to  give  eigenf unctions  expressed  as  a 
finite  number  of  terms  of  the  Hylleraas  series  [19]  •   The  series 
is  defined  in  terms  of  the  coordinates 

s=  tt  +rk      ,   t  =  /r-  £    ,     u  =  r?t  (9.8) 


wh 


ere  YJ     is  an  arbitrary  scaling  parameter.   The  series  has  the 


^0 


advantage  that  wave  functions  composed  of  its  terms  are  automatically 
of  singlet   S   character. 

Hylleraas  [ 19]  found  that  only  six  terms  of  (9.9)  were  required 
to  give  an  excellent  upper  bound.   It  is  these  six  which  we  employ 
in  our  calculations.   We  list  them  here  normalized  to  I6n2  for 
convenience. 


5     24SJ7  b  e  >    h   ~  JT^r?  W  e 


21/17  J  ^  >    r*      *of? 

We  present  first  lower  bounds  determined  by  the  procedure 
introduced  in  Section  VI. 

Elements  Containing  the  Inverse  of  the  Resolvent- 
We  note  from  (9 . 1 )  and  (9.6)  that 

Ei     *    E-2    ,  (9.11) 

so  that  we  may  identify  lower  bounds  to  the  ground  state  by  (6.15). 

Here,  in  order  to  obtain  intermediate  wave  functions  of  the  forn 
(9.9)  we  have  only  to  identify  the  h      of  Section  VI  with  those  of 
(9.10).   Our  lower  bounds  obtained  as  eigenvalues  of  X"'(E~)of 


k] 


orders  one  to  six  require  the  evaluation  of  the  elements, 


each  of  which  is  a  quadratic  in  t-    , 

ty  =  nifty  *  nifty  e"  f  n/jtyfE")*.  0..3) 

Vie   do   not    list   all    twenty-one   of    these   but   rather   give    two    typical 
examples.      We   find   for   /   A,^  the   components, 

M°2    =-fJ]T  \V\l508./3/S  -  (t.5-$)[30/0.9*26     - 

?l  \/szo  -  /6soC/.s-f)r  -/9z\  > 


f  66  g. 5055 


to2  -    63° 

1  *2Z  ~  I  92TJ 

and   for  /   /*^the   components, 


'2* 


M°  = 


Zt        253.992/2 


7/    J  ID36   -  (/  5-^110 


72    - 


(9.15) 


/92 


O  (/-<})]  ~  /803  C/'pt  -7/X[2Z5  -360(1-^)1 


hi 


m,  = 


253.93/21 


( 


^\_/80  0  -  20/6(1. 25-cj) 


)]  -2*5 j  j 


Ml4    = 


7  6a 


2  5  3.  39/Zlff 


In  the  above,  Q     is  the  ratio  2/1?    .   The  procedures  for  computing 
these  elements  are  given  in  Appendix  II. 

We  construct  from  these  quantities  eigenvalue  equations  of  the 
form  (6.11)  of  orders  one  through  six  and  determine  their  eigen- 
values (if  any)  by  finding  the  real  roots  of  the  associated  secular 
equations  (6.12).   The  technique  we  employ  for  finding  the  eigen- 
values and  eigenf unctions  is  an  iterative  one  due  to  Ldwdin   [20]  . 
We  give  a  brief  description  in  Appendix  II. 

We  have  noted  in  Section  VI  that  the  real  roots  of  (6.12)  occur 
in  pairs,  with  an  equation  of  order  /*>    having  from  zero  up  to  •» 
pairs.   In  our  calculations  we  never  find  more  than  one  pair  of  real 
roots.   This  is  despite  the  fact  that,  e.g.,  our  secular  equation  of 
order  six  is  a  polynomial  of  degree  twelve  in  the  energy.   In  Figure 
I  we  give,  as  typical,  the  behavior  of  the  roots  of  a  secular  equation 
of  order  two  based  on  the  elements  /,  and  f.    .      We  see  that,  in  the 
region  centered  around  Tf  ~  2   ,  two  roots  are  found  which  merge  near 
7=  /-75  and  1J  -  2.Z0.      Outside  this  range  no  roots  are  found.   Both 
roots  are  lower  bounds  by  (9.11),  but  we  retain  as  our  lower  bounds 
only  those  given  by  t-  /  . 

It  is  of  some  interest  to  note  that  as  the  order  of  the  cal- 


<+3 


Figure  I.   Roots  of  a  secular  equation  of  order  two. 

2      2 
The  curves  labeled  E.  and  E   give,  at  any  value  ofJf  ,  the  ground  states 

of  different  operators  H2(E2)  and  H2(E2). 


kk 


culations  increases  the  curve   ^  /  flattens  out  rapidly  so  that  at 
/ft  =  6   the  lower  curve  is  qu  i  te  close  to  a  s  trai  ght  1  i  ne,  E  /  —  £->  ~  ~~  '  • 
This  supports  our  discussion  in  Section  VI  regarding  convergence  to 
the  E£   .   The  exact  value  -k    is  attained  at  .£?=*?,  for  any  order, 
since  at  this  point  f,        is  identical  with  H*  i     . 

In  Figure  2  we  give  the  curves  £L  %   obtained  from  our  calculations 
of  orders  one  through  six.     It  is  gratifying  to  observe  the  monotonic 
increase  of  the  lower  bounds  with  order  since  we  have  not  proved  that 
this  must  necessarily  occur  in  this  procedure.   At  points  where  two 
curves  touch  as,  e.g.,  the  point  near  ^  *=  £-0^  for  Jl  j>    and  £l  *      » 
we  find  that  the  higher  order  eigenf unction  reduces  to  the  lower 
order  eigenf unction.   Thus  at  such  points    ^    y  L-  ' 

and  (rC      ( t-      /     have  a  common  eigenf  unct  ion  as  well  as  a  common 
ei  genvalue. 

In  Table  h   we  compare  energies  and  wave  functions  of  our  optimized 
lower  bounds  of  orders  one,  two,  three  and  six  with  those  of  well- 
known  upper  bounds  of  the  same  orders.   We  see  that,  though  individual 
coefficients  and  the  scaling  parameter  differ  somewhat,  corresponding 
wave  functions  have  the  same  general  form. 

For  some  time  two  lower  bound  procedures,  not  connected  with 
intermediate  problems,  have  been  known  and  used  rather  extensively 
in  obtaining  lower  bounds  to  helium  [21,  22,  15].   These  are  due  to 
Temple  [16]  and  Stevenson  and  Crawford  [23]  .   The  procedure  of 
Stevenson  and  Crawford  is  more  difficult  to  apply  than  Temple's  but 

A  calculation  of  order  /w  is  based  on  the  functions  !,,*£,■•  ;/>• 

taken  in  the  order  that  they  appear  in  (9-10). 
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is  somewhat  more  efficient.  '   Actually  it  contains  Temple's  formula 


as  a  special  case  [22]  .   Their  lower  bound  is  given  by 

± 

where  the  quantity  (H$l  frC  €)     -  2  OC  ($OCI  f)]s   minimized  by 
variational  means  and  where  oc  must  satisfy 

oc  -£  £'   +E*  (9.17) 

2 

Optimum  lower  bounds  are  attained  when  the  equality  in  (9.17)  holds. 

In  their  original  paper  [23]  Stevenson  and  Crawford  computed 
lower  bounds  to  helium  using  trial  functions  consisting  of  our 
functions  /,   ;  f,    ,   f2      and  r3     ;  r,       through  / 1      ;  and  a  nine 
term  function  in  which  three  additional  terms  of  the  series  (9-9) 
were  appended  to  the  six  term  function.   We  give  their  results 
in  Table  5.   We  also  compare  our  lower  bounds  with  theirs  and  with 
the  upper  bounds  of  Table  k   by  listing  the  ratios  of  the  corres- 
ponding errors.   We  see  that  our  lower  bounds,  while  quite  inferior 
to  the  upper  bounds,  are  superior  to  those  of  Stevenson  and  Crawford. 
(Actually  our  lower  bound  of  order  six,  -2.909^,  is  very  close  to 
theirs  of  order  nine.) 

Since  we  judge  the  procedure  given  here  to  be  roughly  comparable 
in  difficulty  with  that  of  Stevenson  and  Crawford,  we  consider  the 
data  of  Table  5  to  be  suggestive  of  the  possibility  that  this  pro- 


Kinoshita  [21,  22]  using  a  39  term  trial  function  reduced  the 
error  In  his  lower  bound  by  a  factor  of  two  in  going  from  Temple's 
formula  to  the  procedure  of  Stevenson  and  Crawford. 
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Table  5.   Comparison  of  the  errors  in  lower  and  upper  bounds  to  helium 

Order    Lower  bounds  of    Ratio  of  the  error  Ratio  of  the  error 

Stevenson  and      in  the  lower  bounds  in  our  lower  bounds 

Crawford  (au)      of  Stevenson  and  to  the  errors  in  the 

Crawford  to  the  upper  bounds  of 

errors  in  ours  Table  4. 

(a"s.c./a")  (a"/a+) 


1  -3-5410  3.9  2.9 

3  -2.9481  2.6  13 

6  -2.9215  3.1  12 

9  -2.9089 


Taken  from  [23]  .   The  calculations  of  orders  one,  three  and  si 
could  be  improved  slightly  by  using  a  less  conservative  value  of  a 
[see  (9.17)].   The  ninth  order  calculation  used  a  value  of  a     near 
the  optimum. 
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cedure  may  provide  a  practical  means  of  obtaining  more  accurate  lower 
bounds  than  those  obtained  to  date. 

The  Method  of  Truncations  I 

For  this  procedure  we  construct  the  elements  ^tA-   from  the 
functions  (9.10)  according  to  (7.2),  i.e., 

so  that  our  intermediate  wave  functions  are  expressed  in  terms  of 
the  r*   and  the  /*  as  in  (7-^). 

We  consider  first  truncated  intermediate  Hamiltonians  qC  , 

/h  -  fjZj  ■  •   •  j   6     .      The  eigenvalue  equation  (k.8)    reduces  in  this 
case  to 


(9.19) 

Evaluation  of  the  inner  products  appearing  in  (9.19)  is  elemen- 
tary.  The  formulas  we  use  are  given  in  Appendix  II.   As  an  example 
we  give  the  numerical  results  for  the  equation  (9-19)  of  order  two. 
Denoting  the  elements  in  brackets  by  /V/j  we  find, 
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N      =     875*      +      ***/***'*    +    «>?Z'Z5'/^t'2 
-*-    £>.*  -2.5  -  £<>* 


13 


N/z    =  9.79795927?     +     2  2  &  .  6*  9  Q5/7?*  ! 

-f  -  Eh* 

2  rs     +v'3  (9>20) 

+     3.572I72(>  7f*     ~    Z28.  6  /905/r?*! 

-2.5    -£+* 


/    =  /J>./*5«    /    20  9/666  7/^$      +    *rt2-204.l6667/t¥!'^ 


where  q         is  the  quant  i  ty  (tf+2  '/2  7f  . 

By  multiplying  each  element  /Vxj   of  equation  (3 .$)    by  the 
common  denominator  (T /£"'  y^?5  "^  £ ' J    'we  convert  the  elements 
into  quadratics  in  the  energy    so  that  the  technique  used  to 
solve  (6.11)  may  be  applied  here.   Our  solutions  give,  as  the 
lowest  eigenvalues  of  the  operators  fl     '     )  sn-ZI.'   •   •.  &     » 
the  lower  bounds  in  Table  6.   The  values  of  77       given  are  those 
which  optimize  the  lower  bounds.   The  ei genf unct ion  corresponding 


E  ' 


'j< 


to  C.  is, 


This  operation  introduces  a  spurious  root/",  ~  ~  ¥    into  the 
corresponding  secular  equation  but  does  not  affect  the  already  existing 
roots . 
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4'f     =     <¥*     +   .777°  I     -.3$**  f2     -  .0002   fB 

(9.21) 

+  •1695^    -.O/asfs     +0753  /* 

The  functions  r,    ,  and  the  Tx     are  normalized  to  one  in  (9.21). 

Our  results  in  Table  6  are  not  particularly  good.   In  fact,  our 

best  lower  bound,  £     —  -  3.00  45    is  inferior  to  the  lower  bound 

22 

C.    .      =    ~  3.  OO  O  8      ,  obtained  by  Bazley  and  Fox  [  7  ]  using  this 

procedure.   It  is  clear  upon  comparing  these  two  values  and  noting 
the  slow  improvement  with  Jf      that  in  order  to  improve  these  bounds 
significantly  we  would  have  to  increase  the  index  x       in  our  fC         . 
However,  the  numerical  calculations  increase  in  complexity  very 
rapidly  with  a        ,  so  before  proceeding  to  the  more  involved  cal- 
culations, we  determined  an  upper  bound  to  the  improvement  which  can 
be  obtained  by  increasing  X 

The  upper  bound  depends  on  the  fact  that  the  operator 

a°  =  e:i+;x+ii  +  el[i-  i4>;x+:i] 

=  -4i+:x+;i  -z  [i  -  /+;x+:/] 


(9.22) 


is  larger  than  any  aC         , 

<X'-'>      *     (U°>     ,     (fm    /.*,....).     (9.23) 
Therefore  the  eigenvalues  A  y    of 
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Table  6.   Lower  bounds  to  the  ground  state  of  helium  by  the  method 
of  truncations  i . 


Order  E  '   (au)  Optimum  f 

1  -3.03  2.24 

2  -3.0086  2.05 

3  -3.0084  2.05 

4  -3.0055  2.35 

5  -3.0048  2.40 

6  -3.0045  2.15 

We  quote  this  value  from  Bazley  and  Fox  [7]  . 
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are  upper  bounds  to  the  eigenvalues  of  all  <fC  ,  i.e., 

£"V  *  AZ      ,    C/-/,t.~.  ;'*-*,*....).   0-25) 

The  /ix  different  from  -k   and  -2  are  obtained  from  the  solutions 

of  a  set  of  equations  equivalent  to  (9.19)  except  that  JE %   's  re~ 

F° 

placed  by  L-  O0   . 

We  find  that  the  first  eigenvalue  of  the  ££      corresponding  to 
(rt  is  optimized  for  ff  ~    2.30   and  has  the  value 

a\* =  -  Z.9  6  5  6       .  (9.26) 

Therefore  we  have 

-3.00*5  *  £*',6±  -29656   ,     (9.27) 

On  the  basis  of  (9.27)  we  have  concluded  that  it  would  not  be 
profitable  to  extend  our  calculations  to  larger  values  of  x 

Incidentally,  the  rather  large  difference  between  the  upper 
bound  A  i    and  the  true  ground  state  [L  ,     is  indicative  of  the  well- 
known  fact  that  the  continuum  of  rC     has  a  significant  projection 
on  the  ground  state  wave  function  <J»'/  of  helium  [  \h  ]  . 

The  Method  of  Truncations  II 

As  discussed  in  Section  VII  it  is  not  practical  with  this  pro- 
cedure to  attempt  to  construct  ^4.   which  give  intermediate  eigen- 
functions  expressed  as  linear  combinations  of  the  ^r^  and  the  fa 
of  (9.10).   Our  calculations  with  this  procedure  are  confined  to 
operators  £x    and  are  intended  to  illustrate  the  results  ob- 


5h 


tained  when  the  ^^     are  chosen  so  as  to  lead  to  straightforward 
calculations . 

For  the  operators  oC         the  eigenvalue  equation  (5.^)  reduces 
to 

r  -*-£'••  (9.28) 


-2.5  -  £  ' 


=r  O 


For  our  first  calculation  we  used  the  element, 

A=    uf,     -   0ULJ  «e"7S,  (9.29) 

and  found  as  the  first  root  of  (9.28)  at  7f  =     /.80      the  optimum 
value  -3.8055  for  £  ^    . 


*      -7/5 


p,*    U&      -  ^  *   *  (9.30) 


/,/ 


We  next  tried 

and  found  for  /l  ,   the  poorer  value  -3.898I  at  /7  =■  £.  0 
Finally  we  tried  the  element, 

f.-  uf.   =   ^'|«'ri, 

where 

-7s 


(9.31) 


/„     =     (2*1)  file    '    .  (9.32) 
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Here  we  found  a  lower  bound  of  -3.3^8   at  ff-  '•  '    .  We  consider 
this  value  close  to  the  best  that  can  be  obtained  using  a  single 
element  based  on  simple  functions  such  as  those  used  here.   Further, 
considering  the  poor  results  obtained  with  the  elements  sfi,,      and  jb.^   , 
we  do  not  feel  that  significant  improvement  would  result  in  going  to 
calculations  involving  the  operators  r€        ,  <rC  ,    etc.   Of 

course,  it  may  be  possible  to  construct  elements  which  do  not  lead 
to  unmanageable  calculations  and  which  at  the  same  time  give  good 
lower  bounds,  however,  we  have  found  no  way  of  constructing  such 

elements  for  helium. 

( 

We  give,  as  an  example,  the  e igenf unc t icn  <f> ,       obtained  from 

the  calculation  based  on  the  element  <fl0    , 

~n  r  1  '?lS 

We  see  that  T0      appears  as  the  last  term  in  brackets.   The  wave 
functions  obtained  from  the  other  two  calculations  are  more  in- 
volved  than  (9.33)  and  contain  numerous  terms  in  addition  to  *f>l   and 
the  appropriate  function  T 

We  find  this  procedure  the  most  difficult  to  apply  of  the 
three  used  in  this  section.   The  expressions  for  the  elements  of 
the  eigenvalue  equation  (9.28)  are  even  more  complicated  than 
those  for  the  procedure  of  Section  VI.   To  illustrate  this  we 


In  this  case,  atfj  -  Z  ,    the  method  becomes  identical  to  that 
of  truncations  i.   This  results  from  the  fact  that  atff  -  2  ,   j)}  0 
reduces  to  <fit      and  we  have 

(X~7€'-*)fo  =  />.. 

Using  ^0#  in  a  2f  '"'calculation  Bazley  and  Fox  [7]  obtained  the 
better  lower  bound  -3.29. 
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/rite  out  (9.28)  for  the  calculation  involving  f*, 


■¥ 


■//r'ta  i         i\  «r 


-f- 


256/? 
-4 


4-2147  - 


(/ 5  -})[sr6  '  320(/.5-f)]j  -T(3{/30  -  l2£>(/.5  -f  )J  + 


1*{-7U  +  720(/.5-<})  -^fpy[  35  -/*/ 


V*r 


(9.34) 


+  2 


Ot1   -£*SC/.5-$)]}  +7l[l75  +      jjjfi} 


+ 


60O    +   -zrg 


36750 


tf 


7*  f 


+ 


4/4-20  V2 


256/1*  t"        __     0 
-2.5  -£"" 

In  (9.3^)  Q      is  given  by  Z/*7    and  f      by  (j^t2  V^lf  .      The 

techniques  needed  to  evaluate  the  various  inner  products  in  (9.28) 

are  those  used  in  the  previous  two  procedures  and  are  given  in 
Append  ix  II. 
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AUXILIARY  PROOFS 


APPENDIX  I 

AUXILIARY  PROOFS 

S18 


We  consider  the  vector  space  5  introduced  in  Section  II  character- 
ized by  the  positive  definite  metric  operator  V    and  the  inner  product 
(tylVl^y    >  but  having  its  elements  normalized  according  to 

<4\4y  = ' .  (ai-') 

We  introduce  in  5  a  linearly  independent  set  of  elements  £  ^»,j 
and  consider  the  linear  manifold  //  consisting  of  the  totality  of 
elements  of  5    of  the  form 

We  take  an  arbitrary  element  ^  of  S  and  write  it  as 

4-  4"  +  £ 

(Al-3) 
for  the  purpose  of  showing  that  if  the  Cji,    are  chosen  so  that 


1 8 

This  development,  with  changes  in  notation  to  agree  with  the 

main  text,  is  taken  from  a  seminar  given  by  P.  0.  Ltiwdin  at  the 

University  of  Florida  in  June  of  19&2. 
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then  the  operator  Q      defined  by 

o"4  =  4"  (ai.5) 

is  a  projection  operator  in  D     projecting  on  the  linear  manifold  A*7 
From  (A1.3)  and  (Al .k)   we  have 

/*» 

<Miri4>*E<ftivijn>c;  ,      (A,.6) 

(  f '  -    /,2  ,   •  '  ■  ,  sr>  )  t 
so  that  If  XI   is  the  inverse  of  the  n   order  matrix  with  elements 


From  (Al.7)  we  may  define  (y      as 


0"  =  HZ  I p+)A74<fa\V >       (ai.8) 

<it4»i  ff  " 


since  with    this   definition, 


0  "4  -  E  AL  <mIvi4>* 


(A1.9) 


=  E  c 


£—1       Cx     -^fe-X   4 

The  properties  (2.11)  are  easily  verified  directly  from  (A1.8), 


e.g.  , 

2 


(0")  =2ZJ  /fi>/17,  <t,i  VlfidAu  <filV 

(ALIO! 

=  E  I fii> A?t  <&l IS  =  O"  . 

These,  along  with  the  boundedness  of  yCf    (which  follows  if  V 
is  bounded  in  ft       ),  are  sufficient  to  establish  that  U      is  a  pro- 
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jection  operator  in  S- 

The  inequalities  (2.12)  follow  from  Bessel's  inequality  in  $ 
or  directly  from  the  fact  that  U       is  a  projection  operator  in  S 
and  that  the  subspace  of  is  contains  that  of  */        [25,  p.  68]  . 

We  conclude  this  discussion  by  showing  that  when  the  C j.     are 
given  by  (A1.7)  the  quantity  yO  I V  I  O  /   is  a  minimum.   We  first  take 
O  to  be 

J~  =  ¥    -  L->    A4  fix J  (Al.ll) 

where  the  QA    are  not  restricted  to  be  of  the  form  (Al.7).   For  O 
given  by  (Al .11) 

<HVIi>  =  <4\Vl4>  -g<a.+  fi*IVl4> 

which  from  (Al .6)  can  be  rewritten  as 


(A1.13) 


*****     * 

m 

+  ZJ  Kn+A+IVIaj  fa)  . 

By   adding   and    subtracting     !_,_,  ( CA 4l J V \ C.  A.) 

<61V\6>  =  <4\V\4>  -£Z<C4&IVICjp/> 


we  obtain 


(AI.14) 


which,  since  the  last  term  is  positive  definite,  is  a  minimum  for 


Conditions  Under  Which  the  Transformation 
V'idi0-  €   )         is  Singular 

We  consider  the  equation  (6.2)  defining  the  ^fr*  in  terms  of  the 
linearly  independent  set  (  Lcj  , 

and  seek  the  conditions  under  which  the  transformation  y     ( <r£     ~£  ) 
is  singular.   It  is  sufficient  to  examine  the  transformation  [at    ~~€.J 
since  V       ,  being  positive  definite,  is  never  singular  [5]  . 

We  suppose,  for  convenience,  that  ^rv   's  positive  definite 
and  define  the  functions  Q/   by 

Qj,       ~     (K'-e)^     >     (x*  /,*,-  ;/»).      (A. 2) 

The  necessary  and  sufficient  condition  for  the  4       to  be 
linearly  dependent  is  that  there  exists  a  set  of  coefficients  CLi,    , 
not  all  of  which  are  zero,  such  that 

HaHa      =  0  (A.3) 


*  =  i 


<> 


X°E  •*&    -  e  Ea„l  .  (A.*) 


4.  -I 


1 9     n*o 

If  vC       Is  bounded  from  below  we  make  it  positive  definite  by 
moving  the  zero  of  energy f 

X'-+  7t°+e'     ,    €-*£+■£', 

where    £/>  j £°l 
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Now  neither  side  of  (h.k)    can  vanish  on  account  of  the  positive 
def i ni teness  of  9C      and  the  fact  that  €      is  a  scalar  constant  so 
that  (A. k)    is  just  the  eigenvalue  equation  for  <hC  Therefore,  we 

state  that  the  necessary  and  sufficient  conditions  for  V    (at     ~~  £-  J 
to  be  singular  are  that  some  eigenf unction  r V  of  7€      be  contained 
i  n  the  span  of  f,  t  fg  ,  *  •  •,  Tm   and  that  s  imu  I  taneous  ly  €  —  C-  *    . 
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APPENDIX  II 
CALCULATIONAL  PROCEDURES 
Outline  of  the  Procedures  Used  in  the 
Calculations  of  Section  VIII 

Calculation  of  matrix  elements 

The  elements  of  the  matrices  V     and   1/   of  Section  VIM  are 
given  by 

^/     J  (/y)'  %'  "n  $i)jr2    (/»3   ?<?  ""*   9^^V^(A2.1) 
and 


Z^(A2.2) 


Using    the   relation 

COS  0«     =     n*        r>      ~r>  (A2.3) 

in  equations  (8.10)  and  (8.11)  and  inserting  these  in  (A2.1)  and  (A2.2) 
reduces  the  integrands  to  linear  combinations  of  functions  of  r,  ,  /> 
and  //>  .      Using  a  method  of  Calais  and  LOwd i n  [26]  for  handling  such 
integrals  permits  us  to  write  them  as  linear  combinations  of  the  mul- 
tiple integrals, 

7  4,  +  Ot-/       -X2      «,  az 

J  (x2)  e      Lit  (x2)  Ljlz  (x2)  xTdx 

°    y  (A2.4) 

X  j  <y>)  e    7  LAa<f>LU(i>)Sdf\z*dz  , 

6k 
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where  we  have  let 


X    =  cut  V         ,      y  =  co*rz       ,      Z  =  co*rf2     .    (A2.5) 

Upon  writing  out  the  Laguerre  polynomials  and  carrying  out  the  Z 
integration  these  are  further  reduced  to  sums  of  the  double  in- 
tegrals 


X 

2<t 


C     -x2    2<r+i       C    ~1f 


1 


(A2.6) 


which  may  be  evaluated  in  closed  form, 


Urn,**}*    {Ejr<l>t3H$>y   .  (2f.*,-*-i)  m   (A2  ?) 


(A2.8) 


The  -formula    (A2.7)    is   obtained   from   the  formula    [27   p.    176], 
oe  2  *  L        2 

[e~"X xctx\e'    **„  =  &    ,'.U> 
%  J  /        f  a(a  +6)* 

by  successive  differentiation  with  respect  to  4     and  a. 

A  program  was  written  for  use  on  an  IBM  709  computer  which 
evaluated  the  integrals  (A2.U)  as  linear  combinations  of  (A2.7)- 
These  integrals  were  then  combined  to  give  the  matrix  elements  (A2.1) 
and  (A2.2). 

Inversion  of  V        and  solution  of  the  eigenvalue  equations  (8.15)  and 
(8.17) 

The  ^   were  obtained  from  the  appropriate  submatrices  of  }f       , 
and  the  eigenvalue  equations  (8.15)  and  (8.17)  were  solved  by  computer 
programs  kindly  supplied  by  Mr.  J.  0.  Nordling  of  the  Quantum  Theory 
Project  of  the  University  of  Florida. 


Outline  of  the  Procedures  Used  in  the  Calculations  of  Section  IX 


Evaluation  of  matrix  elements 

The  elements  of  equations  (6.11),  (9-19),  and  (9-28)  require  the 
evaluation  of  one  or  more  of  the  quantities  <K  f,  ft  r    ,  ot        ' 
and  WT    ,  where  f   is  one  of  the  terms  of  (9.9). 

In  [21]  W     is  given  in  the  coordinates  f,     ,  ^   and  f/z    . 
Using  the  relations  (9-8),  the  expression  can  be  transformed  to  the  co- 
ordinates 5,  i     and  u.    with  the  result 


#'«- 


Zt 


s2-t>  Js 


Zst2      J2       +  Jju__  £_      f_2S2t      J*-     (A2.9) 
u(sx-t*)£u<)s       sz  -C*  duds      u(s*-tz)Zudt 

_  ztu    y      ,£_    +  -i-A.   +  -f5^ 

Therefore  &£°   operating  on  f    gives  a  linear  combination  of  terms 
of  the  form 


(s'-tv  s  t     u 


e-'S  , 


(A2.10) 


where  „=   -/  or  0  ,A  i  O     ,/>0     and  ^,  >  -/  .   K  being  u. 
gives  a  single  term  of  the  form  (A2.10)  when  operating  on  /  .   There- 
fore, 7t     ,  being  just  the  sum  of  ?€    and  r     ,    also  gives  linear 
combinations  of  (A2.10).  7C  '"operating  on  T     gives  a  linear  combination 
of  f     and  4>°   •   «^f  (normalized  to  l6tt2  )  is  explicitly 

-zs 
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(A2.ll) 
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With  this  information  it  may  be  seen  from  an  examination  of  (6.11), 
(9.19),  and  (9.28)  that  all  matrix  elements  have  as  integrands  linear 
combinations  of  terms  which  are  at  most  quadratic  in  the  quantities 
(A2.10)  and  in  (A2.ll).   The  elements  themselves  may  therefore  be 
shown  to  be  (see,  e.g.,  [10,  p.  1737])  linear  combinations  of  the 
integral s 

(m/ij.~»)    =  2n*^e~*Ssl<£s 

0 

S  U  (A2.12) 

x$um,<iu^<s2-t'rti'ctt , 

o  o 

where  rn    =-1,0  or  /   ,  4   -  °    .  £  -  ^  .  /y^>   S  -/  and  £      is  ei  ther 

2  ty        or f  +  2  .      These  may  all  be  evaluated  in  closed  form  although 

20 
several  formulas  result.   We  list  them  below. 

fn/i    /         )    -  rrr2(l+2l  +  '»"  +2.)!  !  (A    n) 


WMt£,A*)-2TT       t^+^  +  s        [(2£+,)(2i  +  ~.+ 


2) 

(A2.14) 
/  1 


(2jt  +  3)(2£  +"»+4)  J 
For  m>i  ff, 


(-VA,t.-)***<l:**?;>! 


/m  „  /m-p 


(A2.15) 


For  convenience  in  our  calculations  we  used  functions  normalized 
to   16*    •   To  compensate  we  divided  (A2.12)  by   16«2   in  our  actual 
work. 


-'       a*  =  a 
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,   (2/  +/n,  -2/*X2  J  -2/*  -/) 


(A2.16) 


and  f inal ly 

£-1 

r L 

.jXi      (24     '2y^    '  I) 

A  large  number  of  the  integrals  (A2.12)  were  required  for  the 
various  calculations,  so  a  simple  program  for  evaluating  them  was 
written  for  the  IBM  709,  embodying  the  formulas  (A2.I3)  through  (A2.16). 

Solution  of  the  eigenvalue  equations 

The  elements  of  the  matrix  eigenvalue  equations  (6.11),  (9.19), 
and  (9.28)  are,  or  can  be  made  into,  polynomials  in  the  energy.   We 
therefore  consider  an  n   order  equation 

M  C    =  O       ,  (A2.17) 

where  the  elements  of  A7  are 

4 


"*>  "  f.  "i  U)    •  <+i  -  :>>  ■  ■ :~  >■ 


(A2.18) 


LOwdin  [20]  has  given  an  iterative  procedure  based  on  a  partitioning 
technique  for  solving  such  equations.   The  procedure  was  originally 
developed  to  handle  equations  whose  elements  are  linear  in  the  energy 
but  is  readily  adapted  to  equations  where  the  /*}..   are  more  complicated 
functions  of  £     [  28]  .   We  sketch  here  the  development  of  the  pro- 
cedure when  the  /%y  are   given  by  (A2.18).   We  do  not  consider  compli- 
cations which  arise  In  special  cases  as,  e.g.,  when  there  are  multiple 
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eigenvalues.   These  are  covered  in  detail  in  [20]  . 

We  begin  by  partitioning  /#  and  U     in  (A2.17)  to  obtain 


ribi   rtbt 


=  o  t 


(A2.19) 


where  fi/i     is  one-dimensional  and  fiby\s  0"  -I)  -d  imens  ional  .   Carrying 
out  the  multiplications  we  obtain  the  separate  relations, 


Mn  C,     +    MnCk     =    O     , 

rib.c,  -tMbbCb  =  o  , 


(A2.20) 


which  can  be  recombined  to  give 

( Mn  '  Mlb  Mil  M»,)C,  =  o  .  ^-^ 

If  C,      is  not  zero  we  have  the  relation 

F(E)    =     (M„  -  MhM's's  Mb,)  =    0      ,  (A2.22) 

which  is  equivalent  to  (A2.17),  and  is  suitable  for  solution  by  iterative 
methods . 

In  [20]  it  is  shown,  by  an  application  of  the  Newton-Raphson 
formula  to  Ft C  J ,  that  a  convergent  second-order  iteration  procedure 
for  J-.      can  be  obtained  in  the  form 


E*  -  £  -  UE) 


FU)    ' 


(A2.23) 


where 


F  '(£)  =  J7  F(E)   =   Ml,  -  n(t  Mtb  Mbl 

*-£-  (A2.24) 

+  M,b  Mil  Mbt  M'L  Mt,  - M,b  Mbbrtt, . 
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We  note  that  #7^4  is  easily  evaluated,  having  elements  given  by 


(A2.25) 


The  iterative  process  is  initiated  by  selecting  a  trial  valued 
for  £-     ,  computing  r(£)  and  F (J?)and   determining  a  first  iterate 
*-      ~  *-  x_   is  then  used  in  place  of  C       to  start  a  new  cycle 

which  determines  a  second  iterate  C.      .   The  process  is  continued 
until  the  agreement  between  successive  iterates  is  satisfactory.   The 
eigenvector  belonging  to  the  eigenvalue  thus  found  is,  from  the  second 
of  (A2.20),  given  by 


-7 


C, 


(A2.26) 


evaluated  at  the  eigenvalue. 

A  program,  employing  this  procedure  and  capable  of  handling 

equations  of  orders  up  to  six  whose  elements  are  polynomials  of  degrees 

.   21  .  _ 

up  to  six,   was  prepared  for  use  on  the  IBM  709.   The  program  con- 
tains, as  subroutines,  the  actual  elements  'uV  expressed  as  functions 
of  %       .      Examples  are  the  elements  displayed  in  (9.14)  and  (9.15). 
Input  to  the  program  is  a  set  of  values  of  ff     for  which  solutions  are 
to  be  found  and  a  range  of  energy  which  is  to  be  scanned  for  eigenvalues 


21. 


These  values  are  easily  increased,  being  limited  only  by  loss 
of  significant  figures  and  the  capacity  of  the  machine. 
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